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Abstract 

Among all electromagnetic theories which (a) are derivable from a Lagrangian, (b) satisfy the 
dominant energy condition, and (c) in the weak field limit coincide with classical linear electro- 
magnetics, we identify a certain subclass with the property that the corresponding spherically 
symmetric, asymptotically flat, electrostatic spacetime metric has the mildest possible singular- 
^ ity at its center, namely, a conical singularity on the time axis. The electric field moreover has 

a point defect on the time axis, its total energy is finite, and is equal to the ADM mass of the 
spacetime. By an appropriate scaling of the Lagrangian, one can arrange the total mass and 
total charge of these spacetimes to have any chosen values. For small enough mass-to-charge 
I""! ratio, these spacetimes have no horizons and no trapped null geodesies. We also prove the 

uniqueness of these solutions in the spherically symmetric class, and we conclude by performing 
I a qualitative study of the geodesies and test-charge trajectories of these spacetimes. 

1 Introduction 

In this paper we study static, spherically symmetric solutions of the Einstein-Maxwell system 
(O featuring a single spinless point charg^ Our main results are summarized in a theorem at the end 

of this introduction. The Einstein-Maxwell system of PDEs reads: 

Kf^^-^Rgf,^ = xT^^, (1) 
O dF = 0, (2) 

Tl dM = 0. (3) 

^ It describes the geometry of a spacetime (7W,g) endowed with an electromagnetic field. Here R 

^ is the Ricci curvature tensor and R the scalar curvature of the metric g of a four-dimensional 

Lorentzian manifold Ai. Moreover, F is the Faraday tensor of the electromagnetic field and M 
is the Maxwell tensor corresponding to F(for which we will consider a whole family of choices), 
while T is the electromagnetic energy(density)-momentum (density)-stress tensor associated to F. 
Finally, x = SvrG/c^, with G being Newton's constant of universal gravitation and c the "speed of 
light in vacuum". (In the following we will work with units in which x = 1.) 

Before we discuss the options of choosing the relationship between M and F, we should recall 
that an open domain Ai' in a Lorentzian manifold A4 is (somewhat inappropriately) called static if 
it has a hypersurface-orthogonal Killing field K whose orbits are complete and everywhere timelike 
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'^In a sequel to this paper we will study solutions with spin. 
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in M' . Such a domain possesses a time-function t, i.e. a function defined on it whose gradient is 
timehke, and the vectorfield dual to its gradient is future-directed everywhere. It can always be 
chosen such that = 1. One can use t as a coordinate function on Ai' , and the level sets of t 
provide a foliation of A4' into space-like leaves S^. It follows that the induced metric on the leaves is 
independent of t, i.e. the space St is static. Furthermore, for a "static spacetime" to be spherically 
symmetric means that the rotation group 50(3) acts as a continuous group of isometries on the 
manifold, with orbits that are spacelike 2-spheres, and its action commutes with that of the group 
generated by the timelike Killing field K. For p G A4' , let A(p) be the area of the spherical orbit 
that goes through p, and let r{p) = \J~A{p)~f^ be the radius of a Euclidean sphere with area A{p). 
As long as g(Vr, Vr) > one can use r as a spacelike coordinate function on the manifold, and in 
that case the metric of the spacetime can be put in the form 

g^udx^'dx" = -e^dt^ + e'^dr"^ + r'^{d9'^ + sin^ edcp"^), (4) 

where {9, (p) are spherical coordinates on the orbit spheres, and ^ = ^(r) is a smooth function which 
depends on the choice of relationship between F and M, here called an "aether law" (for short, 
and for historical reasons.) 

The simplest choice of the aether law is Maxwell's M = — * F, in which case the Einstein- 
Maxwell system will be called the Einstein-Maxwell-Maxwell system (EMM) . It is well-known that 
the Reissner-Weyl-Nordstrom spacetim^ (RWN for short) is the uniqu^ spherically symmetric, 
asymptotically flat solution of the EMM system. For the RWN solution, one has 

where qq and tuq are two real parameters. They are in fact integration constants that come from 
solving the radial Liouville-type equation that arises as the reduction of the EMM system under 
the stated symmetry assumptions. Since the spacetime is asymptotically flat, its ADM mass [5J 
is well defined, and it is seen to be equal to the parameter mo. Also, one has a formula for the 
Faraday tensor in the case of RWN: 

F = dA, A = ip{r)dt, (^ = — , 

r 

which suggests, via the divergence theorem, that go is the total charge of the spacetime. 

As discussed in detail in [6], the causal structure of the RWN spacetime depends in a crucial 
way on the ratio \qo\/mQ: When this is less than one, which is referred to as the subextremal 
case, the metric coefficient e^, in addition to being singular at r = 0, has two zeros, namely at 
r± = rriQ ± itIq — . It can be shown that r = r_|_ is the event horizon, the boundary of the past of 
the spacetime's future null infinity, and therefore the spacetime has a nonempty black hole region. 
It is worth noting that the causal structure of the maximal analytic extension of the subextremal 
RWN spacetime is quite rich and complicated, comprising an infinitude of regions, and is plagued 
by the breakdown of determinism due to the presence of Cauchy horizons (cf. [6j.) By contrast, 
the RWN metric in the superextremal regime, corresponding to |go| > w-o has a very simple causal 
structure. The metric coefficient is always positive, (t, r, 9, (j)) is a global coordinate system for 
the manifold, and the only singularity present is the naked one, on the timelike axis r = 0. The 
topology of the manifold is that of minus a line. 

^Generally known as having been discovered independently by Reissner [T] and Nordstrom [2], this spacetime is 
also a member of a whole class of electrovacuum solutions found by Weyl [3] . 
^See Section 6 for a precise statement and proof of uniqueness. 
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In view of the fact that the empirical charge-to-mass ratio of charged particles such as the 
electron and the proton are huge (10^8 and lO^^ 

resp.) many researchers have been tempted by the 
prospect that the superextremal RWN solution is but the simplest example of spacetimes featuring 
one or more point charge^ Such a spacetime we shall refer to as a "charged-particle spacetime." 

One of the key questions that needs to be addressed in this regard is the following: According 
to relativity theory's E = mc^, the proper mass of a charged-particle spacetime should equal its 
energy, which for a static spacetime is expressed as the integral of the time-time component of T 
over a static space-like hypersurface. Is it then possible to attribute some, or all, of the mass of a 
charged-particle spacetime to the energy of the electromagnetic field that permeates itP]? 



the self-energy of a point charge is infinite, and that remains to be the case even when the elec- 
tromagnetic Maxwell-Lorentz field is coupled to gravity via the Einstein equations. This can be 
clearly seen in the case of RWN because the total electrostatic energy carried by a time-slice would 



which is infinitqj unless qq = 0. In general relativity, the principle of equivalence states that it is 
the total energy of a system that interacts gravitationally, i.e. unlike the Newtonian theory, there 
is no distinction between active and passive "gravitational" masses. Thus, the total electrostatic 
energy will always make a contribution to the ADM mass of the spacetime, which in the case of 
the superextremal Reissner-Nordstrom is clearly an infinite contribution. 

On the other hand, the ADM mass is usually interpreted as representing the total energy 
content of the spacetime (recall that in relativity there is no distinction between mass and energy 
(with c = 1)). Ingenious proposals have been made to explain the finiteness of the ADM mass, in 
particular to assume that the point charge possesses a "bare mass" equal to — oo. Such ideas have 
been pursued using renormalization techniques [11] , but these techniques are usually very difficult 
or even impossible to justify rigorousl}|^ 

Another difficulty with taking the superextremal RWN too seriously is the presence of a strong 
naked singularity on the time axis. Such singularities are expected to be "subject to cosmic 
censorship", in the sense that they are believed to be non-generic, and therefore unstable under 
small perturbations. Furthermore, the presence of a strong "eternal" singularity means that such 
a spacetime cannot arise as a solution of a classically-posed initial value problem. For the RWN 
metric the worst part of the singular behavior at r = stems from the contribution of the charge 
qo to the metric coefficient, as can be seen for example from the Kretschmann scalar, which is a 

*See [7], §21.1 for references, and for a catalog of such solutions. 

^Questions about the origin of the mass of charged particles are as old as the theory of electromagnetism itself, 
beginning with O. Heaviside, and continuing with Abraham, Lorentz, Poincare, Mie, Einstein, Fermi, Born, Dirac, 
Wheeler, Feynman, Schwinger, Rohrlich and many others. 

^That is to say, Maxwell's equations with "point-particle-like" sources whose formal law of motion is that of 
Newton's, with the formal force being the Lorentz force. The model is ill-defined unless regularized. 

'^Indeed, this is the same infinity that turns up in the absence of gravity, in flat spacetime, for the self-energy 
of a point charge, and which led Abraham & Lorentz, and later Mie, to look for alternatives to the point-charge 
description. We also note that in the current approach to quantum electrodynamics, the corresponding energy 
integral to the above is still divergent, although less violently [Q, thus even in the absence of gravity, the problem of 
infinite self-energy of point charges is not solved by going over to the quantum description [9].[10]. 

**For example, it is known that a rigorous removal of the regularization using mass renormalization 'Ti^g^j.g — )■ — oo 
is impossible in the case of Lorentz electrodynamics, since the renormalization flow terminates at rn-^^gj.^ = with 
regularization still in place [12j . Similar difficulties exist for the renormalization of RWN [13j . 



One difficulty with such an attribution is that in classical Maxwell-Lorentz electro- 
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curvature invariant equal to the norm of the Riemann tensor: 



Clearly, /C ~ oo as r \ 0. 

To overcome the first of these two divergence problems while retaining the concept of a point 
charge, Max Born [14J proposed to make the Maxwell equations nonlineai]^ This is done by 
choosing a Lagrangian density L for the electromagnetic action S[A] = L((iA) in such a way 
that, in addition to fulfilling the basic requirement of generating a Lorentz- and Weyl- (gauge) 
invariant theory, it coincides in the weak field limit with the Lagrangian of the Maxwell-Maxwell 
systerrp^ 

Lo = --'-FA*F, 

Svr 

while its behavior in the strong field limit is such that it leads to finite total energies. One example 
of such a Lagrangian density is the well-known one-parameter family proposed by Born and Infeld 

L/3 = - Vl - /34 * (F A *F) - /38(*(F A F))2], 

for /3 > (in the notation of |18j). which even leads to finite limits of the field strengths at the 
location of the charge. 

What is perhaps less common knowledge -though it should be equally well-known- is that a 
nonlinear aether law also has the power to reduce the strength of the spacetime singularity that is 



present when the electromagnetic field is coupled to gravitj This phenomenon was first noticed 



by Hoffmann [20] , who initially claimed that a solution of the Einstein-Maxwell-Born- Infeld system 



that he had found, was free of all singularities In the years since the publication of Hoffmann's 



paper there have been several attempts at finding static electrovacuum spacetimes which are free of 
all singularities l2ll[22l[23l|2l[23[2ni|271[28], either in the fields or in the metric, while at the same 
time various "no-go" results have been announced [291 EH] , that seem to show that such solutions 
cannot exist. One of the goals of this paper is to take steps towards dispelling the confusion that 
seems to persist about this subject. 

Our approach here is to characterize aether laws that not only feature finite self-energies for 



point charges, but also lead to electro statu^^ spacetimes with the mildest form of singularity pos- 
siblep^ Thus we will initially allow the aether law to have the most general form possible, and let 
the above requirements, of finiteness of the self-energy and mildness of the singularity, as well as 
other criteria which may arise during the course of analyzing the corresponding solutions, dictate 
the final form that it should have. In particular we will prove the following: 



^Born picked up on the program initiated by Mie [15], who however did not want point charges. See |18l 119] for 
an exceUent account of the development of these ideas, and for their author's key contribution to this program, 
which seems to have been discovered by Schwarzschild [16j . 
^^Even though, as we will prove in this paper, it cannot quite eliminate the singularity. 

^■^Hoffmann's enthusiasm for nonlinear electrodynamics was not dampened even after it was pointed out to him 
by Einstein and Rosen (See [21| . fn. 15), that a mild singularity remains at the center of his spacetime. 
^^See Section 3 for a precise definition. 

^*In fact, we will see that the remaining mild singularity in the field and the metric is of the point-defect type. 
There are indications that this kind of singularity may be "just right" for the Hamilton-Jacobi equations to provide 
a law of motion for those defects |18l I19| . In this way it may turn out then that the "underwater stone" of nonlinear 
electrodynamics (in the words of [S^) is a gem after all! 
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THEOREM. (Informal version) For any aether law which (a) agrees with that of Maxwell in 
the weak field limit, (b) is derivable from a Lagrangian, with an energy tensor that satisfies the 
Dominant Energy Condition, and for which (c) the corresponding Hamiltonian satisfies certain 
growth conditions (to be made precise below), the following hold: 

• There exists a unique electrostatic, spherically symmetric, asymptotically flat solution of the 
Einstein-Maxwell system {^-(^ with that aether law, the maximal analytical extension of 
which is homeomorphic to M'^ minus a line. It has a conical singularity on the time axil , 
which is the mildest possible singularity for any spherically symmetric electrostatic spacetime 
whose aether law satisfies (a) and (b) . No other singularities are present in the spacetime. 

• A generalization of Birkhoff's Theorem shows that this solution is unique in the spherically 
symmetric class. 

• The electrostatic potential is finite on the axis of symmetry, which can be identified with 
the world line of a point charge. The electric field has a point defect at the location of the 
charge. The field has finite total electrostatic energy, which is equal to the ADM mass of the 
spacetime. The mass of the point charge is thus entirely of electromagnetic origin, i.e. it has 
no bare mass. 

• By scaling the Hamiltonian appropriately, one can arrange the total mass and the total charge 
of this solution to have any chosen values. The mass-to-charge ratio of the spacetime enters 
as a natural small parameter, measuring the departure from the Minkowski spacetime. 

• For small enough mass-to-charge ratio, there are no horizons of any kind and no trapped null 
geodesies in this spacetime. 

• The analysis of geodesies and test-charge trajectories shows that the naked singularity at the 
center of this spacetime is gravitationally attractive (unlike the case of superextremal RWN). 

The rest of this paper is organized as follows: In section 2 we introduce the Lagrangian for- 
mulation of electromagnetics. Section 3 gives the equations satisfied by an electrostatic solution 
of ([T])-([3]), with an arbitrary aether law. In Section 4 we assume spherical symmetry as well, and 
obtain the general solution to the equations. Section 5 is devoted to the study of the singularities 
of this solution. In Section 6 we state and prove a Birkhoff-type uniqueness result for these space- 
times. In Section 7 we give the precise version of our main result, and in Section 8 we carry out 
the qualitative analysis of geodesies and test-charge orbits. 

2 Nonlinear Electrodynamics 

Let (A^,g) be a 4-dimensional Lorentzian manifold. Let denote the bundle of p-forms on 

M. By an electromagnetic Lagrangian (density) we mean a mapping Lem defined on sections of 
the vector bundle A^(-^) A^(-^) which takes its values in (see [38] for details.) Thus 

if a is a 1-form on Ai and f a 2-form, then Lem(a, f) is a 4-form on Ai. The electromagnetic action 
is by definition 




^^That is to say, the limit as the radius goes to zero of the ratio of the circumference to the radius of a small 
spacelike circle going around the axis, exists but is not equal to 2tt. 



5 



where P is a domain in . A critical point of S with respect to variations of a that are compactly 
supported in V is called an electromagnetic potential A in V, 



6S_ 



= 0, 

a=A 



and the exterior derivative of it is the electromagnetic Faraday tensor F = dA. The Maxwell tensor 
M is by definition 

M 



(5) 

a=A,f=F 



in the sense of evaluation, i.e. the object on the right, when evaluated on a variation f G 
^p(A^(-^)) = A^(^p-^) a four- form is equal to M A f . The source- free Maxwell equations 
are the Euler-Lagrange equations for stationary points of the electromagnetic action S, and are 
equivalent to the system 

dF = 0, dM = 0. 

It can be shown |38j that the only Lorentz-invariant gauge-invariant source-free electromagnetic 
Lagrangians possible are those of the form 

Le^(a,f) = -^(x(f),y(f))e[g], (6) 

where e[g] is the volume form on Ai induced by the metric g, £ is a real- valued function of two 
variables, and x and y are the electromagnetic invariants 

x(f) := -1 * (f A *f) = y(f) := 1 * (f A f) = * 



Here * is the Hodge star operator * : /\^{Ai) — )• /\^ ^{A4) with respect to the g metric, defined by 



*C/xi.../X4_fc — , ,^ ^ ''£[g]!^l...!^fcMl...At4-fc- 



1 

Note that for fc-forms on a Lorentzian 4-manifold, ** = (—1)^"^^. 

Furthermore, conservation of parity implies that i{x, y) = i{x, —y), so that if we assume that i 
is a function of its variables, then 

£y{x,0) = 0. 

By an aether law we simply mean a particular function i{x, y) as the Lagrangian density func- 
tion, which determines the way the electromagnetic vacuum interacts with the spacetime geome- 
try For example, conventional (linear) electromagnetics corresponds to the choice I = x made 
first by Maxwell. 

Using * one defines a dot product on /c-forms, by 



cr • r = - * (cr A *r) = ^<T^'-'''=r^i...^fe . 



We also set \a\^ := a ■ a even though this is not necessarily positive. It follows from ^ that 
*Lem = £ and 

*M= — = 4F + £y*F. 



^''Traditionally, aether law referred to the relationship between tensors M and F, similar to the constitutive 
relations of elastodynamics relating stresses to strains for the medium. In case of a Lagrangian theory, this is given 
by Isl, and thus the choice is that of a particular Lagrangian. 
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The energy tensor T corresponding to the Lagrangian density function £ is a symmetric 2- 
covariant tensor field on A4 defined by 

di 

which in the case of an electromagnetic Lagrangian yields 

T^, = 2 (42F^aF,^ + iy^Y^y^ *F,^) - ^^J = *M,^ - g^,£. (7) 

Recall that if both of the following hold, the energy tensor T is said to satisfy the Dominant 
Energy Condition: 

• T^pY^y^^ > for every future-directed timelike vector Y, 

• The vector — Tt'Y'^ is future directed causal when Y is a future-directed causal vector. 

The first of these two is called the Weak Energy Condition. There is also a Strong 
Energy Condition: 

'T^,-^g^,trT) Y'^Y''>0, 

for all future-directed timelike Y. One can prove [39j that the Dominant Energy Condition is 
satisfied for this field theory if and only if 

4 > 0, I- x4 - y^y > 0. 

Note that under the above assumption, tr T = —4{i—xix—yiy) < 0. We also note that it is possible 
for the Strong Energy Condition to be violated in nonlinear electrodynamics (even though, as we 
shall see, it can hold for nonlinear electrostatics.) 

Next we define the following one-forms, which provide a useful decomposition of the Faraday 



and Maxwell tensors. Let K be an arbitrar}!^ non-null vectorfield on Ai, and let^^ 



iKF (8) 

iK *F (9) 

iK *M = i^iB + £y^ (10) 

iK**M = -iKM = -£y(B + i^^, (11) 



where denotes the interior product with the vectorfield K, e.g. {iKF)u = K^F^,y. 

For the Maxwell Lagrangian £ = x one has S) = 2; and ^ = S). A general aether law will specify 
D and 55 as functions of and S), or the other way around. 

Let 

X:=g(K,K). 



^'^In the next section we will assume that K is a timelike Killing field for g), but the definitions in this section 
are independent of that. 

'^^Note that strictly speaking, only if K were the unit tangent vectorfield to a timelike curve in M (the world- 
line of an observer) would we be justified in calling these the (flattened) electric field, magnetic induction, electric 
displacement, and magnetic field, respectively. 
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Thus X > wherever K is spacehke, X = where K is null, and X < where K is timelike. 
From the decomposition of F in terms of <B and 55 it follows that 



We also obtain that, by ([T]) 



T{K,K) = (B-^-Xe, (13) 



while from ([T]) it follows that 

R(K, K) = (T(K, K) - f trT) = (*B • + Xl). (14) 

For K a timelike vectorfield, we can define two electromagnetic Hamiltonians (partial Legendre 
transforms of the Lagrangian with respect to either G; or 5S~ 
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H{ia,S)) := \X\-'m ■ — - £ = \X\-^^ ■ Sj - i, 

H(s,^) := -\x\'^^- — + e = \x\-'^(s-T> + e. 

In the definition of H (resp. H) we need to think of 53 (resp. C;) as given by the aether law. Also, 
the factors of X in these will disappear in the next Section, once the inner products are re-expressed 
in terms of a different metric on space- like slices which is conformal to the one induced by g. 
We then have 

T(K,K) = |X|^(5S,D), R(K,K) = \X\H{(B,S)). 

The Weak, Dominant, and Strong Energy Conditions, respectively, take the following form flO] in 
terms of the Hamiltonians 

H>0, H>\H\, H>0. (15) 

It is easy to see [JT] that H is actually only a function of the three scalar invariants that one 
can form out of ^ and using the metric. More precisely, 

where 



and that h G C^(M^) is such that its gradient lies on the future unit hyperboloid in M3, i.e. 

hi - hi- hi = 1, h, > 0. 
Similar results hold for H, in terms of the invariants ;U, tu, a defined analogously in terms of S, OS: 

H{^,D) = h{fj,,vu,a), 

where 

2^(D-D + 5S-5S), w := • 53, a := ^(D • D - *B • 55), 
and likewise h £ C^{]R.^) with gradient lying on the future unit hyperboloid. 



^These can also be defined for a spacelike vectorfield, but there will be some sign changes. 
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3 Electrostatic Spacetimes 



Let K be a timelike hypersurface-orthogonal Killing field for the spacetime , g) . Let X = 
g(K,K) < and define 

e = dX. 

Let 

denote the metric induced on the quotient Q of under the symmetry generated by K, and let 

7 = |X|h, 

so that 7 is also a Riemannian metric on Q, conformal to h. Since K is assumed to be twist-free, 
the quotient Q can be identified with a spacelike hypersurface in M.. Let (j;*), i = 1,2,3 be an 
arbitrary coordinate system on Q. Then the line element of g is given by 

ds'^ = Xdf + \X\-^-fijdx'dxK 

g 7 1 

It follows in particular that V -V = X V ^^'^ ^'^y vectorfield V such that CyiV = 0. The 

Einstein-Maxwell system in the static case reduces [IQ] to the following set of equations: 



V-(^e) = =iH{^,^) (16) 

7 



V-C^S) = (17) 

V-(i53) = (18) 

RicMi,- - \l^jRb] = 2^e,e,- + ^(£^1), + *B,io,-) - 7,, (4^|e|^ + ^i/((S,io)) , (19) 

where to close the system one has to remember that 



with b being the operation of lowering indices with respect to the 7 metric. 



The equation ( 19 ) confirms that we can interpret the above as a system of Einstein equations 
for the manifold Q coupled to the fields (e, We further recall that Maxwell's equations for 

F and M also imply that (2; and S) are exact 1-forms 

^ = dip, S) = dif). 

Thus the above is a system of equations for the three potentials X,ip,ip, and the metric 7 of the 
quotient manifold. 

We also observe that the parity conservation assumption about the Lagrangian ^y(x,0) = 
implies that, away from singularities, !l) = whenever G; = 0, and likewise = whenever *B = 0. 
This fact, together with the invariance of the equations under interchanging T) with 55 and C; with 
io, implies that the system of equations in the magnetostatic case £ = is formally the same as 
that in the electrostatic case = 0, even though (it turns out) the restrictions on the Hamiltonian 
under which meaningful solutions can be obtained are differenlPj 



"For examples of regular magnetostatic (magnetic monopole) solutions, see [26] . 
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In this paper we are confining ourselves to the electrostatic case, where u} = 'nj = 0, t = i/ and 
a = 11. If we define 

r]{ij) = h{u, 0, u), C(/") = H^J'^ 0, n), 

then in the electrostatic case ^ = §f" = ^ ~ f§ ~ C'(m)®- We can moreover express 

rj in terms of This is because in the electrostatic case *B = ^ = 0, and thus 

+ = 0. £), = (!-. If), = 

SO that we have 

77(z.) = 2K'(^)-C(m). (21) 



In terms of the reduced Hamiltonian ^ the Dominant Energy Condition (15) takes the following 
simple form 

C'(/i)>0, C-K'(/^)>0, V/i>0, (22) 
while the Strong Energy Condition reads 

2K'(/^) - C(/^) > 0, V^>0. (23) 

We note that in the electrostatic case, u = —x and y = 0, thus ?7(z^) = —i{x, 0) = —£{—h', 0). On 
the other hand the function ^ can be obtained from i via a Legendre transform: Given £ = i{x, y) 
let f{t) := —i{—^t'^,0) and let /*(s) = sup^[st — /(t)] be the Legendre transform of /. Then it 
is easy to see that Cil^) = f*{\/^)- As an example, here are the Lagrangian density function 
originally proposed by Born \T^, and its two reduced Hamiltonians: 

eB{x) = Vl + 2^-l, 7]b{i^) = 1-Vi^, Csif^) = ^1 + 2^-1. 

Let 

^=log(-X). 

With u = ^\d(p\'^^ and S = r]'(iy)dip, the electrostatic Einstein-Maxwell system then becomes 

V -d^ = 2e-^7]{u) (24) 
V -(e-^D) = (25) 



(26) 



4 Spherical Symmetry 



If we assume that the spacetime (and the electromagnetic field) in addition to being static is 
also spherically symmetric, there will be a further reduction in the Einstein-Maxwell system. In 
particular, using the area-radius coordinate r we may rewrite the metric 7 as follows 



-fijdx'dx^ = dr^ + e^r'^{d9'^ + 



sm 



where now ^ = ^(r), and also (p = (p{r), D = D{r)dr. From (25), 



=V -(e-^D) = -^dMr^e-^D) = -^dr{r''D), 
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and thus one obtains that 

^ c , 

D = -^dr, 

where c is an arbitrary constant. On the other hand, since 



^='l^l7 = ^- (27) 



From (24) we now obtain 



1 d ( f 



r^e^ dr \ dr 
Change variable to n = ^ to obtain 



e^r — I = 2e ri{v). 



We can now use (21), integrate (28) twice, change the order if integration and recompute the 
kernel to obtain a formula for ^ function of r: 

,ii'^^j'^i^iriM%, (29) 

r r Jj. fi r'^ 

where c,cf,cf' are arbitrary constants, and = The requirement that the solution be asymp- 
totically flat now forces c" = 1. On the other hand, setting C = should give the Schwarzschild 
solution, hence c' = — 2mo where mo is the mass parameter in the Schwarzschild metric. 

We can also find an expression for the electrostatic potential (/? as a function of r. Recall that 

dr 

^'{r)dr = d^ = <a = C'(/x)D = cC'(/x)^, 

and thus 

f°° , dr' 
<^(r) = cj C'(^)^- 

Since in the Maxwell-Maxwell case Cif^) = comparison with the RWN solution shows that c = qo, 



the charge parameter in RWN. Finally, a direct computation shows j40j that (26) is identically 
satisfied. 

Thus we have the following generalization of RWN to nonlinear gravitoelectrostatics: 

r r 

('(^^)^■^ ^ = 1^.- (31) 

These simple and elegant formulae seem to have first appeared in [23], with a different deriva- 
tion, although special cases of it were known long before (see below.) Note that only the reduced 
Hamiltonian ( makes an appearance in them. The prospect of generating exact solutions to the 
Einstein-Maxwell system with interesting and desirable properties, just by inserting a judiciously 
chosen into these formulae has proven to be irresistible to many theoreticians. In particular we 
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should mention the solution found by Hoffmann [22] to the Einstein-Maxwell- Born-Infeld system, 
which corresponds to the following Hamiltonian: 



Cb/(/^) = + (32) 



Another early example is the solution found by Infeld and Hoffmann |21j . where they made the 
following choice 

Cm(/x) = log(l + ^), (33) 

and obtained a metric which was completely smooth and free of all singularities! Their work was 
followed up by Rao [22j, who attempted to find a large family of actions leading to such solutions. 

Infeld and Hoffmann however may also be the first of many investigators who have made the 
mistake of picking a Hamiltonian that is not admissible because it cannot arise from a Lagrangian: 
One has to remember that C(/f^) is the electrostatic reduction of the Hamiltonian H{^,Ti), which 
in turn is subject to the following restrictions: 

1. It must in the weak field limit agree with Maxwell's choice for the aether law. 

2. It must correspond to an energy tensor that satisfies the dominant energy condition. 

3. It must be the Legendre-Fenchel transform in C; of a Lagrangian density i = ^), i.e., it 
must be convex in S). 

It follows from the above that the function ^ is subject to the following requirements: 
(Rl) hm^^o ^ = 1. 

(R2) C'ip) > and C - K'(Ai) > for all > 0. 
(R3) C'ip) + 2/uC"(/i) > for ah /i > 0. 



The Hamiltonian ( |33[ ) proposed by Infeld and Hoffmann violates the third condition above. There- 
fore it cannot arise from a single-valued Lagrangian. 



The condition (R3) is equivalent to insisting that for the above solution (31), u = —■^{'■p'ir))'^ 
be monotone decreasing in r (note that /x = qQ/(2r^) is always monotone, independent of the choice 
of a Lagrangian.) It was shown by Bronnikov et al [29j that this monotonicity requirement of u 
rules out the possibility of having electrostatic solutions with a regular center (i.e. no curvature 
blowup at r = 0) that are at the same time Maxwellian in the weak field limit. The same argument 



applies to show that the three conditions above are incompatible with having a regular centei^^ 



Therefore spacetimes corresponding to ( 30||3T ) must have some kind of a singularity at r = 0. 



In view of this fact, fantastic claims about existence of singularity- free point-charge metrics in 
nonlinear electrodynamics, which every now and then appear in the literature, should be viewed 
with a healthy dose of skepticism, and the "Hamiltonian" involved should be examined carefully, 
for it may violate one or more of the above conditions. 

We note also that (R3) implies in particular that the Strong Energy Condition, which we men- 
tioned before can be violated in nonlinear electrodynamics, nevertheless holds in the electrostatic 
case, since < 2^C" + (' = (2/iC' — (")', which upon integration on [0,fj] and using (Rl) implies 



that (23) holds. This will have important consequences, as we will see below. 



^^Another possibility is of course, not to have a center at all [26;. The topology of such a spacetime however, does 
not seem to lend itself to the point-charge concept. 
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Interestingly, a question that should have been addressed long ago, but wasn't, is this: Since 
the above requirements rule out solutions that are everywhere regular, what is then the mildest 
singular behavior allowed by them? In the next section we characterize those static spherically 
symmetric point-charge metrics that have the mildest form of singularity possible at their center. 



5 Singularities 

5.1 Singularities in the metric 

We begin by calculating the curvature tensor of a static spherically symmetric metric of the form 

ds^ = fdt^ - f-^dr^ + r^dn^, 
with / = f{r). The nonzero components of the Riemann tensor are 

1^0101 = //" + /'^) R-0202 = R-OSOS = — R-1212 = — R-lSia = ^ff, R-2323 = ^(1 — /^)- 
The indices here refer to the rigid frame {u>^^^} defined as follows: 

a;(°) = fdt, J^'^ = f~^dr, J"^^ = rdd, = rsin^#. 

Thus the Kretschman scalar is, with X = — 

IC^ = X"^ + r-^X'^ + r-\l + Xf. (34) 



It is evident from ( 34 ) that for there to be no spacetime curvature singularity at r = it is necessary 



and sufficient that 1 + X{r) = 02{r ) as r — t- More generally, the Kretschman scalar /C will 



blow up like r ° if and only if 1 + X(r) = 02{r^ 

For spherically symmetric, electrostatic spacetimes, 



2 

X = —1 H — m(r) 
r 



where 



m(r) ■=1710- C Jj4^ ^''^d^' 



is the mass function. 

Thus we see that /C will blow up at least like if m(0) ^ 0, as it is for example in the case 
of the "negative mass" Schwarzschild solution, where mo < and qo = 0. For the superextremal 
RWN solution, the situation is much worse since m(0) = — oo. Since our goal here is to characterize 
solutions which are as mildly singular as possible at the location of the charge, we may start by 
requiring m(0) = 0. Now, since 

"^(0) = - y-"\{y)dy = mo - kol'/'/c 

and mo is an integration constant, it is always possible to meet this requirement as long as 

/•oo 

k ■■= 2~''/' / y'''^\{y)dy < oo. (35) 
Jo 

From now on we will add this to the list of requirements that the reduced Hamiltonian must 
satisfy: 



■^■^For a continuous function / : — >■ E, integer k and a > we say that / = Ofc(r°) if lim, — s-o r-' ""d? f /dr-' exists 
and is finite for j = Q, . . . ,k. 
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(R4) 

roo 

pr"^ /'^C{lJL)dll < OO. 



Note that this new restriction implies the fohowing: If C(m) is assumed to grow hke a power /x", 
then we must have a < 3/4, which of course rules out the Maxwellian case. 

Having made such a choice of the integration constant mo, we now observe that 

rriADM = m{oo) = |go|^^^^c- 

This means that the mass of this spacetime is entirely of electrical nature. Moreover, by an ap- 
propriate scaling of the aether law, namely ip{x,y) := /3~^^(/3*^x, /3*^7/), which for the reduced 
Hamiltonian amounts to using 

= /3-'C(/3V) 

in place of it is possible to "fit" the ADM mass and total charge of this solution to the empirical 
mass and charge of any particle. This is because J,^^ = so for a given pair of numbers (q, m) 

setting go = Q and 

(36) 

M 

will result in m(oo) = M, while at the same time the scaled version of an admissible Hamiltonian 
function C, will remain admissible. In this connection it is worth mentioning that if one carries out 
this procedure for the Born-Infeld Lagrangian, with mo and set to the mass and charge of the 
electron, the value for the scaling parameter /3 thus obtained coincides with the value originally 
proposed by Born [18] . 

Once mo is chosen as above, the mass function can be rewritten as follows: 



"^(^) = / r'Hr' = M- Tcpi^A r''dr'. (37) 



f 

2r'^J ' \2r 
We can now obtain some estimates for m(r). Recall the second part of (R2): 

C(/^) > 

Integrating this inequality on [0, fi] and using (Rl) we easily obtain 

C(/u) < fi V/i > 0. 
This gives the following lower bound for the mass function: 

m(r) > M , 

which in turn gives the following bound on the metric coefficient e^: 

< 1 + 



These two estimates are clearly only useful for large r. In fact, assuming slightly more than (Rl), 
one can turn these into large-r asymptotics for m and e^. Namely, let us assume that 
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my 



Substituting into the second expression for the mass function in (37) we obtain 
m(r) = M--^ + ( ^] , e« = l + -^ + 0' 



As advertised, these asymptotics are identical to those of the RWN solution. 

We now need to establish the small-r behavior of the mass function. Clearly, m{r) = 02{r^) 
if and only if C, goes to a constant as — )• oo. On the other hand, we recall (R3) on i.e. the 
convexity requirement. It is equivalent to 

{^'/'C{^)y > 0. (38) 

Integrating this on an interval [^O) fA ^ fixed fiQ > 0, we obtain 

C(Ai) > Ci + C2^, 



with Ci, C2 7^ constants depending on ^q. Therefore the reduced Hamiltonian must grow at least 
like fj,^^^ in order to satisfy this requirement. In particular, it cannot go to a constant at infinity, 
hence there will be curvature blowup at r = no matter what aether law is chosen, as anticipated 
in He]. 

Assuming now that the reduced Hamiltonian grows at the slowest possible rate, namely like 
/x^/^, we obtain that there must be a conical singularity at r = where IC blows up like r~^, and 
that there are no horizons for this metric. In order to do this rigorously, we need to make the 
growth condition more precise: 

(R5) There exists positive constants J(^,K(^,L(^ depending only on the profile (" such that 



2^1/2 ^ -^y^^ - 2^1/2- 

We should here pause to mention that an example of a reduced Hamiltonian satisfying all five 
requirements (Rl-5) is the one of Born-Infeld (32). Many other examples can easily be constructed 
by considering smooth, monotone increasing, concave functions of /i that behave like /i for small // 
and like Cy/Ji for large /i. 

Assuming (R5), from (37) we have. 



m{r)<-e'r, A:=—^, e:=j^. (39) 



This right away implies that there will be no horizons as long as e is small enough, since 

-X = l-^-^>l-Ae'>0. 

r 

Combining (R5) with our previously obtained bounds for the reduced Hamiltonian, 

max{0, - K(;} < ({fi) < min{/i, J^^}, 

max{0,^-^}< C'(^) <min{l,^}. 
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We can now use this to obtain a lower bound for the mass function that does not degenerate near 
r = 0, namely 

, , Ae"^ Kc . ,,,, 
m{r) > —r-r - ir^r^ , (41) 



as well as the following small r asymptotics for m: 

m(r) = r :rr +02(r^), B^ := — ^. 

^ ^ 2 2m^ ^ 3/| 

Consequently, 

-X = e^ = + 6^2 + 02(r2) 

with a = \/l — Ae^ and 6 = ^e^- This in turn implies that there is a conical singularity at r = 0, 
because the line element of the spacetime metric is Xdt^ — X'^dr'^ + r'^{d9'^ + sin^ OdcfP'), and the 
coefficient of c?r^ at r = is ^, which is greater than one. In fact, introducing standard Cartesian 
space coordinates (x^,x^,x^) near r = 0, with r = we see that the line element in 

these coordinates is 



xdf + 



1 , \ x'^x^ 



dx^dx^ , 



\X\ J r2 

thus the metric has no continuous extension at r = unless a = 1. We therefore need to take the 
axis r = out of the spacetime manifold Ai, which gives it the topology of minus a line. 

We also note that, given a profile C) the deficit angle of the conical singularity is proportional 
to e = m/|q|. This will be quite small precisely when the empirical charge-to- mass ratio of the 
particle to which this solution is being fitted is large, which happens to be the case for the electron 
and the proton, etc. This means that in the study of these metrics it is permissible to consider the 
small e regime. 

The last observation we would like to make about the metric before moving on to discussing 
the electric field, is that the metric coefficient = —X is monotone increasing. The easiest way to 



see this is from (28). Recall that the assumption (R3) implies that the Strong Energy Condition 
is satisfied, and thus ry(i^) > 0, so that is convex as a function of u = K We have already 
established that < < 1, e^\u=o = 1, lim^-^-oo = < 1. Thus cannot have a local 
maximum or a local minimum at a finite n, and must therefore be monotone decreasing in n, hence 
monotone increasing in r. We will see that this fact, which is in great contrast to the behavior of 
the same metric coefficient in the RWN spacetimes, has important consequences, in particular for 
the trajectories of test particles. 

5.2 Singularities in the electric field 

In the spherically symmetric, electrostatic case, the only nonzero component of the Faraday tensor 
is Frt = '■p'{t) and the one- form = ii^Y = dip. We have 

One easily computes that 

3 A 

sgn(Q)^(0) = -e, sgn(Q)99'(0) = -^e^- 
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Moreover, using (40), we have 



r IqI Ae^] , ^ „ ^ f Ce^ 2 

max< — >< sgn(q ) ip Ir) <min<0,---; 1 ^-r 

[ 2m J \ j-r \ J y 2m M'^ 

with C := 2L^//^. In particular, 99 is monotone decreasing, and is asymptotic to the Coulomb 
potential as r — )■ c«. Moreover, (^'(0) 7^ 0, and thus ^ becomes undefined at r = 0. More precisely, 
since <B = ip'(r)dr and dr is a unit covector whose direction is undefined at r = 0, it follows that 
the covectorfield <B is of finite magnitude and undefined direction, i.e. has a point defect at r = 0. 
<t is otherwise smooth. 

We now check that the total electrostatic energy is finite, and in fact is equal to the ADM mass: 
By virtue of ([T]), the energy tensor T is divergence free, i.e. V^T^j^ = 0. Let := — T^j^K^ where 
K = 9f is the timelike Killing field of (7W,g). It follows that 6P = *d * P = and thus by the 
divergence theorem *P is a conserved current, i.e. /^^ *f* independent of t. With coordinates 
{t, r, 6, (p) as before, K = (1, 0, 0, 0) and we calculate that the only nonzero component of *P is 

(*P)i23 = -g°°Po\/-detg = ^T(K, K)V-detg = H{^, S)^/-detg = C{f,y sin 9. 

Thus, defining the total electromagnetic energy carried by the slice St to be 

E:= / *P, 

47r 



St 



we see that for the particle-spacetimes under study 



1 pZiT pn poo poo 

^ = ir / / (*P)i23 drdOdcP = / Citiydr = M, 
4vr Jo Jq Jo Jq 

as promised. 



6 Uniqueness in the Spherically Symmetric Class 

Birkhoff 's celebrated theorem [3IJ states that any spherically symmetric solution of vacuum Ein- 



23 



This was 



stein's equations Ti^,^ = is locally isometric to a region in Schwarzschild spacetime 
generalized to the electrovacuum case by Hoffmann [36J, thereby proving the uniqueness of the 
RWN solution in the spherically symmetric class. Many other extensions and generalizations have 
followed since, see |37j and references therein. Here we prove that the charged-particle spacetime 
whose existence we have established in the above, enjoys the same uniqueness property. 

To begin, we recall the definition of spherical symmetry: A Lorentzian manifold (A^, g) is locally 
spherically symmetric if every point in Ai has an open neighborhood V such that 



■^^This result is often paraphrased inaccurately as "any spherically symmetric vacuum solution of Einstein's equa- 
tions must be static," which makes it no longer true in general [33] ■ The theorem was discovered first by J. T. Jebsen 
|32| . whose proof appeared two years before Birkhoff 's. Jebsen's proof however, contains an error [33]. There is a 
more general result, by J. Eiesland ,34^, on necessary and sufficient conditions for the existence of an extra Killing field 
for spherically symmetric (not necessarily vacuum) spacetimes, a preliminary version of which was also announced 
|35| before Birkhoff's book but the final paper appeared two years after it. Birkhoff's theorem is a corollary of 
Eiesland's result. 
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1. There is a linearly independent set of vector fields {Uj}^^^ on V which generate a faithful 
representation of the Lie algebra of the rotation group SO{3), and whose orbits are two- 
dimensional and spacelike. (Note that the orbits do not need to be contained in V). 

2. Cij.g = in V for i = 1, 2, 3, i.e. Uj are Killing vector fields for the metric. 
Next we recall the following (see [M] for a proof): 

Proposition 6.1. The necessary and sufficient condition for local spherical symmetry of a smooth 
Lorentzian manifold (A^,g) is that every point in M has a neighborhood in which there exists a 
system of local coordinates {t, r, 6, (p) such that the line element of the metric in those coordinates 
has the form 

ds^ = -A^{t, r)dt^ + B^{t, r)dr^ + C^{t, r){de^ + sin^ (42) 
with A, B, C smooth functions of {t, r] , 

In this case, the three Killing fields generating spherical symmetry are 

Ui = sin (pdo + cot cos (pd^p, U2 = cos (pdo — cot 9 sin (pd^, U3 = d^f,, 

and we can check that these vectorfields satisfy the commutation relations of the Lie algebra of the 
rotation group: 

[U,,U,]=6f^.Ufc, 

where Cijk is the completely antisymmetric 3-symbol. We also note that [Ui, Z] = if Z is either 
dt or dr- Moreover, we can express 80 and in terms of (Ui,U2): 

(9e = sin 0Ui + cos (/)U2 , 5^ = tan0(cos(/)Ui — sin(/)U2). (43) 

Suppose now that F = Y ^ydx^ f\dx'^ is a spherically symmetric 2-form defined in a neighborhood 
V oi a, spherically symmetric manifold (A^,g), i.e. Cjj^F = for i = 1,2,3. Let X, Y,Z be any 
three vectorfields on A4. For any 2-form F we have the identity 

(£xF)(Y, Z) = XF(Y, Z) - F([X, Y], Z) - F(Y, [X, Z]). 

Thus, for a spherically symmetric F, we obtain the following 

1. UiF(at,a,) = 0, i = 1,2,3. 

2. UiF(Uj,Z) = ef^.F(Ufc,Z) for i,j = 1,2,3 and Z G {dt,dr]. 

3. UiF(Uj-,Ufc) = for distinct, and U,F(Ui,Uj) = e^F(Ui,Ufc) for i = 1,2,3. 
Thus, from item (1), deF{dt, dr) = d^F{dt, dr) = which implies that 

Fidt,dr) = fit,r). (44) 

From item (2) we have UiF(U3,Z) = -F(U2,Z), U2F(U3,Z) = F(Ui,Z) and U3F(U3,Z) = 0, 
giving a differential equation for F(U3,Z), and solving that we obtain that for each choice of Z 
there is a function a{t, r) such that 

F{d^,Z) = a{t,r) sine. (45) 

■^"^Jebsen's error was that he had assumed C{t,r) = r. This can only be achieved if VC is spaceHke, thus for a 
complete proof one had to consider three other cases as well, of VC being timelike, null, or zero. 
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On the other hand, item (2) also imphes U3F(Ui,Z) = F(U2,Z), U3F(U2,Z) = -F(Ui,Z) 



which gives a differential system for F(Uj,Z), for i = 1,2. Solving it, making use of (43)and (45) 
one obtains that for each choice of Z there is a function b = b{t, r) such that 

F{de,Z) = b{t,r). (46) 

Finally, in a similar manner item (3) gives rise to a differential system for F(U3,Uj), i = 1,2, 
which upon solving yields that there is a function c(i, r) such that 

F{de,d^) = c{t,r) sine. 

We have thus shown that a spherically symmetric 2-form F on a spherically symmetric manifold, 
with coordinates (t, r, 9, (j)) as above, must have the following form 

Ftr = f{t,r) Fte = bi{t,r) Fr,e = b2{t,r) 

F t(j) = ai{t,r) sin Fr^ = a2{t,r) sin0 F g^j, = c{t,r) sin 9. 

Now, assume that F is a closed 2-form: dF = 0. Using the identity CxF = ixdF + dixF valid 
for any vectorfield X and any tensor F we obtain that for a spherically symmetric closed 2-Form 
F, 

diu,F = 0, i = 1,2,3. 



Writing the resulting differential equations for the components of F as found in ( 47 ) one then sees 
that the only solution is 

ai = 02 = 6i = &2 = 0, c{t, r) = c, 

with c an arbitrary constant. Thus we have proven the following Lemma (stated without proof in 
fc36i): 

Lemma 6.2. On a spherically symmetric manifold M with coordinates {t, r, 9, (j)) the only nonzero 
components of a spherically symmetric Faraday tensor F are 

Ftr = f{t,r), Fe4> = csin9, (48) 

where f is an arbitrary function and c an arbitrary constant. 

Note that the constant c here is the total magnetic charge: 

where S is any spacelike 2-sphere in that contains the origin. 

We can now compute the energy tensor T. Recall that the line element of the metric g has the 
form (42). The nonzero components of the dual tensor *F are therefore 

AB — 
*Ftr = ^c, *Fg^ = ^^f{t,r)sin9, 

and the dual to the Maxwell tensor *M = i^F -\- iy * F is computed to have components 

AB / \ 

* Mtr = 4/ + cly-^, *Me^ = f c4 - ^J^f) sin^. (49) 

We now quote the main theorem in [3l] (slightly reworded to match our notation): 
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THEOREM 6.3. (Eiesland, 1925) The necessary and sufficient conditions that a locally spheri- 
cally symmetric Lorentzian manifold with line element 



ds^ = -A^dt^ + B'^dr'^ + C^{d9'^ + sin^ 

A, B and C being arbitrary functions oft and r, and C not a constant, shall admit a one-parameter 
group of isometrics generated by a vectorfield K = kQ{t,r)dt + ki{t,r)dr are 

A^Gi = ^'drCdtC, (50) 

ab{g;-gI) = ^'c[A\drCf + B\dtCf], (5i) 

where G'^ = IV^ — ^RS'j^ is the Einstein tensor and ^' = ^{C) is an arbitrary function of C, or a 
constant. If the above conditions are satisfied, then the Killing field K is, up to a constant multiple, 
given by 

We are now in a position to state our uniqueness result: 

THEOREM 6.4. Let {M, g, F) be a locally spherically symmetric solution of the Einstein-Maxwell 
system ([Ip-Q). There exists a system of local coordinates {t,r,6,(f)) on M. in which the line element 
of the metric g takes the form (42). If^C is spacelike in A4, then there exists a Killing vectorfield 
K which is hypersurf ace- orthogonal and timelike everywhere in M, i.e. the solution is static as 
well, and in the case of no magnetic charge it is thus isomorphic to a region in the electrostatic 
charged-particle spacetime found in our Section 4- 



PROOF: The existence of the coordinate system is guaranteed by Prop. 6.1 By Lemma [6.2| the 
Faraday tensor F in these coordinates has the form (48). On the other hand, by virtue of ([T]) we 
must have G'^ = T'j^. The relevant components of the energy tensor T are easily computed from 
(|48j) and ^ to be: 

= T^ = ^2^2 (^^^f + '^y~(y2^ ~ ^' '^i = 0- 

We thus see that the solution satisfies the conditions of Eiesland's theorem, with ^ a constant. We 
can take the components of the Killing field to be /cq = —drC/{AB) and ki = dtC/{AB). In that 
case 

g(K, K) = -{drC/Bf + {dtC/Af = -g-\VC, VC), 

and thus K is timelike if and only if VC is spacelike. Since dK is easily seen to be proportional to 
dt A dr, the twist of K vanishes: K A dK = 0, and therefore K is hypersurface-orthogonal, i.e. the 
solution is static . It would then be a matter of changing the coordinates (t, r) into new coordinates 
(t', r') that satisfy 

dt' Bt' Br' Br' 

Kt' = ko^ + k,^ = l, Kr' = ko^ + k,^=0, (52) 
Bt or Bt Br 

which is solvable by quadratures, in order for the line element of g in the new coordinates to have 
the form 

ds"^ = -P^dt'^ + Q^dr'^ + R^{d9^ + sin^ 
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where P,Q,R are functions of r' only, and R is not a constant. At this point we may again use 
that Gq — = 0, where the indices now refer to the new coordinates = t', = r' . Prom the 
definition of G^ we compute (see e.g. [M]) that 

G[] = (Q3 + 2Q'RR' - QR!^ - 2QRR") , G\ = ^^2]^ {PQ^ - PR'^ - 2RR'P') . 

The equaUty of the two is now easily seen to imply that there exists a constant A such that 
PQ = XR'. Thus letting r = At' and taking {T,R,9,(f)) as the new coordinates, the metric takes 
the form 

p2 \2 

ds^ = -^dT^ + -^dR^ + R^(d6^ + sin^ edcp'^), 
A^ -P^ 

which is the form we assumed in Section 4 in order to derive the static spherically symmetric 
solution we found. 

Remark 6.5. The above uniqueness result can most likely be strengthened by considering the 
remaining cases, of VC being timelike, null, or zero, finding in each case what the solution reduces 
to, similar to the treatment of the cosmological vacuum solutions in [37J. We do not pursue this 
here, however. 

7 Precise Statement of the Main Result 

We are now in a position to give the precise version of the main result: 

THEOREM 7.1. Let ( : M^" — ^ IR+ be any function satisfying the following conditions: 

1. ((^) = ^ + 0(^5/4) as n^O. 

2. C'(/u) > and C - K' > for all fi>0. 

3. C'(/u) + 2^C"(/^) > for all /i > 0. 

4. /c:=/o°^/.-n(^)d/.<oo. 



5. There exists positive constants ixq,J(^,K(^,L(^ such that Jc_^/Jx — < ({fi) < J(^^/Jj., and 
1^ - ^ < C'(/^) < ^jjri hold for fj. > fj.Q. 

Let M > and Q ^ be two fixed real numbers, and let e := m/|q|. Define 

1 . o3/2/. 
CM^) := ^C(/3V), /3 := 

Let i be any electromagnetic Lagrangian density function with the property that the electrostatic 
reduction of its Hamiltonian is the function C,p . Then the following hold: 

1. The Einstein-Maxwell system ([^]{^ with the electromagnetic Lagrangian density h = —*£ has 
a unique electrostatic, spherically symmetric, asymptotically flat solution, the maximal ana- 
lytic extension of which is a Lorentzian manifold (A^,g), called a charged-particle spacetime. 
Ai is topologically equivalent to minus a line. There exists a global coordinate system 
{t,r,0,4>) on Ai with the property that K := dt is an everywhere-timelike Killing field, r is 
the area-radius of rotation group orbits, r = is invariant under K, and {9, (p) are standard 
spherical coordinates on S^. 
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2. Any region in a spherically symmetric electrovacuum spacetime where the area-radius of the 
rotation group orbits has a spacelike gradient, is necessarily static as well, and in the case of 
no magnetic charge is isometric to a region in the corresponding charged-particle spacetime 
(A^,g) defined above. 

3. The line element of the metric g has the form 

ds^ = -e^dt^ + e-^dr^ + r^{de^ + sin^ Odcf?), 

where 

m(r) is the mass function of the spacetime Ai. In particular, m{0) = 0, m is increasing in 
r, and the ADM mass of M is m(oo) = M. The mass function moreover has the following 
asymptotics: 

, , Ae^ Bh^ 3 ^ , , Q' ^ n A 

m[r) = ——r — - — + 02[r ) as r — )■ U, m[r) = M — - — as r ^ oo, 



2 2m2 ^ ' 2r \r^^ 

where A := ^J'' o.'^d B^ := Consequently, grr\r=o = (1 ~ Ae"^)^^ > 1 and thus the 

metric has a conical singularity at r = 0. 

4- The mass function satisfies the following bounds 

Ae^ Bh^ 3 Q^l ^ , , ^ . f 



max s r ^-r , M > < mlr) < min < r, M 

2 2m2 ' 2r j - ^ ^ - 12' 

Consequently, for e small enough, gtt < everywhere, which implies that there are no horizons 
in A4, hence the conical singularity at r = is naked. 

5. r = is a curvature singularity, where the Kretschman scalar blows up like r^"^ . Any spher- 
ically symmetric electrostatic spacetime homeomorphic to minus a line, whose reduced 
Hamiltonian satisfies the first three conditions listed above for will necessarily have a sin- 
gularity at r = which is at least this strong. 

6. The spacetime Ai is endowed with an electromagnetic field F = dA, where A = —ip{r)dt, 

roo / „2 



P fQ^\ ds 



is the electrostatic potential, and E := e ^^"^dip is the (flattened) electric fielc , The potential 
if is smooth, monotone decreasing, and has the following asymptotics 



sgn(Q)ip(r) = r + 0(r ) as r — )• 0, ip(r) = ^ as r ^ oo. 

2 2m r \r-^ J 

Thus the total charge of the spacetime is Q. The electric field E is smooth everywhere except 
at r = where it has a point-defect, i.e. its magnitude has a finite limit but its direction is 
undefined. 



I.e., E — ij^F where K is the unit vectorfield in the direction of K. 
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7. The total electrostatic energy carried by a time slice is finite and is equal to M, therefore the 
mass of the spacetime is entirely of electric origin. 

8. Radial null geodesies fall into the singularity and are thus incomplete in one direction, while 
null geodesies with nonzero angular momentum are infinitely extendible in both directions. 
For e small, these are deflected by the singularity by an amount that is proportional to e^ , 
and there are no trapped null geodesies. 

9. Similarly, timelike geodesies and test-charge trajectories can only reach the singularity if they 
are radial, while non-radial ones have either bound orbits or escape orbits. The singularity at 
r = 0, in contrast to the one in the superextremal RWN solution, is gravitationally attractive. 

The only items in the above Theorem that we have not yet proved are the last two, regarding 
geodesies. These will be established below: 



8 Geodesies and Test-Charge Trajectories 

Consider the following Lagrangian density, defined on a velocity bundle (see |38j for definitions): 

^(q, v) = ^ga/3(q)v"v^ + eAA(q)v\ 

where A is an electromagnetic vector potential defined on A^, q G Ai and v G TqTW. The 
corresponding action 

^[q] = / L(q,q)ds 

is defined on curves q(s) = (q"(s)) in Ai, where the dot represents differentiation with respect to 
an affine parameter s, which is related to arclength parametrization ("proper time" for timelike 
curves) r by r = ms. The stationary points of this action are, for e = and = 1,0,-1 
respectively timelike, null, and spacelike geodesies of the spacetime. Moreover, stationary points of 



the action with e ^ and = 1 represent possible worldlines of a "test charge"!^ of mass m and 
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charge e in the spacetime (A^,g) permeated by the electromagnetic field F. The Euler-Lagrange 
equations for these geodesies are 

where D/dr is covariant differentiation operation on tangent vectors, and F is the electromagnetic 
field tensor. 
Let 

dL 

Pa = = ga/3(q)v^ + eAa(q) 



^®This is a fairly standard treatment of test charges, cf. [431 144| . [45) . [46) and [47[ . Notice however that it is not 
without conceptual problems: If one thinks, e.g. as in [46] and [47], that the charged-particle spacetime in question 
is the spacetime of an electron, then one cannot possibly "test" an electron with an electron, i.e. a particle which has 
charge comparable to that of the charged-particle spacetime can by no stretch of imagination be considered a test 
particle, since its effect on the geometry of spacetime cannot be ignored, while on the other hand if the spacetime is 
to represent an elementary particle then of course no particle of smaller charge exists to test it with! The situation 
considered here is thus merely a cartoon, and proper treatment of this subject is postponed to a future paper, where 
we plan to consider charged-particle spacetimes featuring two point charges. 
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be the canonical momenta, and % the corresponding Hamiltonian to L defined by 



-H = p„v" -L = ^g"^(q)(p„ - eA,(q))(p^ - cA^(q)) = ^|p|2 - eA • p + '^\A\\ 

The geodesic equations in Hamiltonian form are q = ^ and p = — The first constant 
of motion is H itself, and from the normalization condition discussed above we have that along 
solutions H. = —]^m^. Moreover, J^p, q) is a constant of motion iff {J^'H} = 0. 

Let us now take g to be a point-charge metric, with empirical charge Q and mass M (that is 
entirely of electric origin), g is the spherically symmetric electrostatic solution to the Einstein- 
Maxwell system with a nonlinear aether law, corresponding to a reduced Hamiltonian C,p that is 
subject to the restrictions (Rl-5) discussed in the previous sections but otherwise arbitrary, and 
let F be the corresponding electromagnetic field permeating this spacetime. The line element of g 
thus has the form 

dsl = -e^dt^ + e-^dr^ + r^{d9^ + sin^ edcf)'^), 

where {9,(j)) are spherical coordinates on the standard unit sphere S^. Moreover F = dA where 
A = (p(r)dt, and the functions ^(r), (p{r) are given in terms of C by 

« = 1-?!!^; rcJ^)s'ds-, (53) 







r 

^ = qI Cpi^J^-^ Qp{x) = p-\{p^x), (54) 



with /3 as in (36). 



With coordinates (x") = {t,r,0,(f)) thus chosen, we have 

Pt = -e^i + tip; pr = e~^r; pg = r^9; p^ = sir? 6<j), 

and 

^ = J \-e-Hpt - t^f + ehl + \pI + 2 ^ 2 a ] ■ (55) 
2 [ r^sm y J 

Since dH/dt = dT-l/dcp = we get two more constants of motion, pt and p^, and we call the values 
that they take on a solution the energy E and angular momentum about the cj) axis 

Pt = E; p^ = 

A fourth constant is also found, by observing that 

d 2h d dU 
—r 9 = —pg = -— = ——3- cost*. 
ds ds oO r^ sm^ 6 

Thus, if one can arrange it so that 9 = Tr/2 when 9 = 0, then 9 = as well, so that 9 will stay 
constant at 7r/2. But since we still have the freedom of choosing the axes for the two spherical 
coordinates, we can always choose the (p axis in such a way that the plane 9 = 7r/2 is the plane 
through the origin that contains the geodesic's initial velocity vector, so that ^ = initially. We 
have thus shown that all geodesies of these particle-spacetimes are planar, and for any single 
geodesic we can always assume that it is contained in the equatorial plane 9 = Tr/2. 
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Now, evaluating (55) along a solution we obtain that (— e ^{E — cip)"^ + e ^r^ + m^)r^ = — <I>^, 



which can be used to find an equation for r: 



$2 

2 



r 

With four constants of motion the equations are reduced to a first-order system of ODEs. We 
proceed to study this system, first for timelike and null geodesies (e = 0) and then for test-charges 
(e/0). 

8.1 Geodesies 

For timelike (resp. null) geodesies, c = and m = 1 (resp. m = 0). The equations are 



i = -e-^E, r = ±\lE^- {^+mAe^, ^=\^- (56) 



We observe that for $ = the geodesic is radial, with r = ±V E"^ — xn^e^. Thus, since < e« < 1, 
for E < ma there is no solution. For ma < E < m the geodesic does not have enough energy to 
escape, and will fall radially towards the singularity r = where curvature blows up. The geodesic 
reaches the singularity at a finite parameter value, since the function under the square root has 
only a simple zero. This is also the case for geodesies with E > m (and therefore for all null 
geodesies), the only difference being that they can reach infinity, and in fact have a well-defined 
asymptotic velocity: (dr / dt)\r=oo = \/l — xn^/E^. All radial geodesies are therefore inextendible 
in one direction and thus incomplete. 

Let $ > now. It is more convenient to reformulate the equations in terms of a reciprocal 
radial variable, and to measure length in units of M. Let 



M 

X := —. 

r 



From the r equation in (56) one easily obtains 

deb ±1 



^52_(2;2+^2)(l_2^M(f))' 

where 

Mm uE M 

^^=^' ^■=^' ''-=1^1' 
and M{y) is the normalized mass function 

^(y) ■= ^nmr / c(/")M-'/'d/^ = 1 - .^YYTir / c(M)/u-'/'d/^, 



which, in view of (|40|), satisfies the following bounds 

max{l-iy,|-^}< M{y) < min{l, A} 

max{-i,2^}< M'iy) < min{0, ^ + |^}, 



1 -At ^ n,r//..\ ^ ^;^rn -A , SS^ ^ 
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where 

A.-^>1, B 

are constants that depend only on the profile (. Note that for the RWN metric, MjiwNiv) = 1 — ^U- 
Let 

f^Ax) := {x' + j') (l-2a;M(J)). 

We note that for null geodesies, 7 = 0. In that case, it is not hard to see that fo,e{^) = x^e^^^^^^ is 
monotone increasing, for e small enough, thus there is a unique xq > such that 5 = /o,e(a;o), and 
we must have x < xq along any null geodesic. Hence ro := m/xq is the perihelion for the geodesic, 
i.e. the closest it can get to the singularity. Furthermore, since fo,t has no critical points, it follows 
that there are no bounded orbits for null geodesies, so that all null geodesies with $ > can be 
extended in both directions to an infinite value for the parameter, and the events corresponding to 
those infinite values are points at infinity. This can be easily seen from the f equation 

as = — , = tor ro < r < oo, 

^JE^ - f e« 

where s is any affine parameter along the null geodesic. Since under the square-root has only a 
simple zero at r = ro and goes to a constant as r ^ oo we see that the integral of the right-hand- 
side of the above is finite over any finite subinterval of [ro, oo); and that this integral diverges only 
when the upper limit is infinite. 

We thus have {dcp/dx)'^ = l/(/o,e(xo) — /o,e(a;)). For the Minkowski space, fo,ei^) = 
xq = S, thus the solution is ^(x) — 0(0) = sin'~^(x/(5), which as expected is a straight line at 
a distance 1/5 from the origin. Since the particle-spacetime is always asymptotically Euclidean, 
$/(m£') = 1/5 is the impact parameter, the distance from the central singularity of the initial 
asymptote of the geodesic. Geodesies with 1/(5 = go straight to the central singularity, while 
geodesies with 1/(5 > will be deflected by it, as we will see below: 

We find it more convenient to parameterize the family of null geodesies in the equatorial plane 
by their reciprocal perihelion xq. Let (pg denote the corresponding angle to xq. We then have 



XQ 



±dx' 



The integral is convergent at the upper limit since as we said above denominator has only a simple 
zero there. Letting (^± denote the two values for cf) obtained at x = 0, we have 

r" (l:r 
(f)+-(P_ = 2 === =: 7r + d(xo). 

^0 ^/foA^o) - JoA^) 

The quantity d{xo) represents the total deflexion of null geodesies due to spacetime curvature. 
For the Minkowski space, one easily computes that d(xo) = 0,Vxo > 0. By contrast for a static, 
spherically symmetric spacetime whose line element is —a dt^ + {l/a^)dr^ + r^idO"^ + sin^ ed(j)^) 
with < 1, so that it has a conical singularity at r = 0, we have /(x) = a^x^ and as a result 
d{xo) = Tr{l - 1) > 0. One can thus say that the null geodesies are "bent" by the gravitational 
attraction of the conical singularity. Below we will establish the lower bound d{xo) > ce^ when 
xo is sufficiently large, or equivalently, when ro/M is small. For large ro/M on the other hand, the 
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geodesic stays far away from the singularity and thus its qualitative behavior is the same as in the 
RWN spacetime, as analyzed in [45j and more extensively in [4Sj §40]. 

As to the lower bound, we have, using that c?x^ < /o < c?x^ + 6^, with = 1 — Ae^ and 
62 = 2B\^ 

1 u \\ 1 _i axQ 1 _i 1 

{7r + d{xQ))> / —===== = -sm > - sm 



2 



^fixo) - a?x^ a ^/(^ a ^\ + (hjax^f 



Let < c < A'njA: be fixed and assume x^) > ^^^_^ e. We then have 



/4- 



and thus 




which is the desired lower bound. 

Next we consider time- like geodesies with nonzero angular momentum. One sees that /^^^(O) = 
7^ > 0, f!y^^{0) = —7^ < and that the bounds for M{y) obtained above imply that f-y^eix) grows 
like for large x. Thus once again geodesies with (5 > 7 (i.e. those with £^ > m) will have 
a perihelion and are extendible to infinity in either direction. In contrast to the null geodesies, 
however, there are bounded orbits for time-like geodesies. This is because /^^^ will have at least 
one critical point. Let /=„ := min^. f^^^{x) = f^^^[x^). It then follows that there will be no solution 
with 5 < ^/Jl, and that for ^/Jl < (5 < 7 the equation /^^^(x) = 5"^ has at least two solutions, 
corresponding to the perihelion and the aphelion of a bound orbit. There will also be at least one 
stable circular orbit, for 5 = \/J*. 



8.2 Trajectories of test charges 

For the RWN metric, a detailed study of test-charge trajectories was done in [IT]. We are not aware 
of a comparable study done for the generalized RWN metrics with a nonlinear aether law, other 
than the particular case of the Born-Infeld law [49] in the black hole regime, and some preliminary 
discussion of the general case in [50]. Here we are going to analyze the qualitative behavior of all 
the test-charge trajectories in a given particle-spacetime with a small enough mass-to-charge ratio 
e. 

Recall that the orbit of a single test particle of mass m and charge e, in an electrostatic, 
spherically symmetric particle-spacetime of mass M and charge Q satisfies the following system: 



In addition, the orbit lies in a plane, which is taken to be the plane 6 = 7r/2. The system is clearly 
invariant under the simultaneous sign change of -E, e, ^ and the independent variable s. Therefore, 
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it is enough to consider the case E >0. Let x := M/r as before. It then follows that 

where e = |^ as before and we have introduced new parameters 

e/m E ^ $ 

p:=——, k:=—, a := , 

Q/M m mM 

and where N{y) is the normalized electrostatic potential: 



Jo 



which in view of (40) satisfies the bound 



, 5 A] ^ ( 3 A C 



Moreover, N'{y) = C'{I^y^/2), so that 

max |o, ^ - ^1 < N'{y) < min |l, 



A 

y2 



Note that for the RWN metric, NRWNiy) = U- 

K and A are proportional to the energy and angular momentum of the test charge. The ratio p 
is positive in the case where the particle-spacetime and the test particle have charges of the same 
sign, and negative if they have the opposite sign. Setting p = will reproduce the results obtained 
above for timelike geodesies. Our goal is to obtain a bifurcation diagram in the k, p, A parameter 
space for the above system, for a fixed small value of e. It is clear that the diagram will be invariant 
under the simultaneous sign change of k and p, and thus it is enough to restrict our attention to 
K > 0. 

Let 

gp,^,e{x) :=K-pN , hx,e ■■= (1 + X^x^) (l - 2xM (- 

Consider first the case of orbits with zero angular momentum: <I> = A = 0. The motion of the test 
charge is then radial. We recall that ho^eix) = e^^^^^^ Thus, /io,e is convex, decreasing, /io,e(0) = 1, 
/iQg(O) = —2 and lima;_i.oo ho^^{x) = = 1 — Ae'^. 

On the other hand, gp,ft,e{0) = k, g'p^i^^^{0) = -p/e^ and that Yim^^oo 9p,n,e{x) = k - 3p/2. 
Accordingly, there are three main parameter regimes to consider: 

Case 1. p < 0. The two charges Q and e thus have the opposite sign. The function gp^K,e 
is increasing, and has a horizontal asymptote. Thus (7p ^ e ™ust have the same properties, and in 
particular it will be concave and increasing, while we have already established that /io,e is convex and 
decreasing. It then follows that there will be no trajectory with /t < a + 3p/2, for a + 3/9/2 < k < 1 
the trajectory is bounded, and corresponds to a test charge that cannot escape and falls radially 
inward and into the singularity, and for k > 1 the test charge has enough energy to escape to 
infinity. 

Case 2. p > 0, k > 3p/2. The two charges thus have the same sign, and the test charge has 
relatively high energy. The function gp,K,t is decreasing, has a positive asymptotic value, and is 
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convex. Thus gj, wih have these same three properties as well. In order to find the number of 
intersections of the graphs of ffp^^.e ^^^^ ^o,e; one needs to compare the values of the derivatives of 
these two functions at x = and near x = oo. we have 

while, using the large y asymptotics established above for the two functions M{y) and N{y), we 
have 

25p,«;,e(x)5p^«^,(x) - h'o^.ix) = -pe^ '^^'^ (;^) ^ ^ as x ^ oo. 

Thus we need to distinguish two subcases: 

(2a) K < I p = (Note that this requires p^ < 2e^ /Z, or equivalently, e/m < a/273, 

which does not allow the test charge to have a charge-to-mass ratio typical of actual elementary 
particles.) In this case there will be a critical value of k, say G (3p/2,a -|- 3p/2), such that the 
graphs of the two functions gp^Kc,e ^0,e are tangent to each other at a positive x = Xc- Thus for 
K < Kc there will be no trajectories, for k = Kc there will be a stable equilibrium, corresponding 
to the test charge remaining at rest with respect to the singularity, for Hc < k < a + 3/5/2 the test 
charge shuttles back and forth between a perihelion and an aphelion, which grow further apart as 
K is increased, for a + 3p/2 < k < 1 the perihelion coincides with the singularity, and for k> 1 the 
aphelion is at r = oo, i.e. the test charge is allowed to escape. 

(2b) K > e'^ / p: The two above-mentioned graphs can have no point of tangency, irrespective of 
K. Moreover, there is no intersection at a positive x, and hence no trajectory, with k < 1 (i.e. for 
E < m). For 1 < k < a -|- 3p/2 there is a single intersection, which corresponds to the perihelion 
of a hyperbolic trajectory that avoids the singularity, and for k > a + 3p/2 the test charge in one 
direction falls into the singularity and in the other direction is allowed to escape to infinity. 

Case 3. yO > and k < 3p/2. The charges have the same sign and the test charge has relatively 
low energy. The function gp^K,e has a negative asymptotic value, thus the graph of (7p^K,e ^i^^ be 
tangent to the x-axis at some point x = x*, which will be where the global minimum is achieved, 
and is asymptotic to (k — 3p/2)'^. There are no trajectories with 3p/2 — a < k < 1. For k < 
min{l, 3p/2 — a} the trajectories correspond to falling-in test particles. For each 1 < k < 3p/2 — a 
there is an avoiding trajectory as well as one that falls in. For 3p/2 — a < k < 3p/2 all trajectories 
avoid the singularity. 

Based on the above analysis, the n-p parameter plane is divided into several regions as in 
Figure [!} The shaded region in this figure is a forbidden zone: there are no trajectories for (k, p) 
in this region. Region I corresponds to shuttling orbits, region II to test particles that are trapped 
and fall into the singularity. Particles with parameters in region III can also fall into the center, 
but have enough energy so that in the opposite direction they are allowed to escape to infinity. 
Region IV corresponds to hyperbolic trajectories that avoid the singularity altogether, while each 
point in region V corresponds to two trajectories, one that falls into the center and another that 
avoids it. 

Next, we consider charged trajectories with nonzero angular momentum, i.e. A > 0. The only 
thing that changes in the above analysis is the behavior of /iA,e) which now grows like x^ for large x, 
and has a global minimum at a positive x. As a result, there will be no trajectories falling into the 
singularity. The trajectories in fact will be similar to classical Keplerian trajectories, divided into 
hound and escape (scattering) orbits. The following diagram (Fig. [2]) shows the various parameter 
regimes for each type of orbit. 
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K 
1 




Figure 1: Bifurcation diagram for trajectories with no angular momentum 



K 
1 




Figure 2: Bifurcation diagram for trajectories with positive angular momentum 

The shaded region in Fig. [2] is a forbidden one. The parameters in regions labeled E give rise to 
bound orbits, while those in regions labeled H give rise to escape (scattering) ones. One observes 
that, in sharp contrast to the case of superextremal RWN studied in |47j . there is no evidence here 
of gravity being repulsive in the vicinity of the naked conical singularity at the center of these 
particle-spacetimes. On the contrary, the singularity appears to be attractive. 
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